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摘 XB: 蒸 散 (Evapotranspiration,ET) 是 生态 系统 水 循环 中 的 重要 一 环 ,决定 了 生态 系统 水 分 和 热 


量 传输 。 从 区 域 义 度 对 蒸 散 及 
环境 因素 对 其 的 影响 机 制 , 有 助 于 合理 利用 、 分 
环 的 影响 提供 参考 。 基 于 生 : 
型 ,验证 模型 在 研究 


t 2k i (Transpiration , T) fu 2& /£ (Evaporation, ,已 ) 组 分 进行 量化 ,认识 
配水 资 
态 系统 生产 力 模 拟 (Boreal ecosystem productivity simulator, BEPS ) 模 
区 域 的 适用 性 ,量化 1981 一 2018 年 内 蒙古 半 干 旱 区 的 ET 及 其 组 分 的 变化 情 


源 ,为 研究 气候 变化 对 区 域 生 态 系统 水 文 特 


况 , 并 对 其 进行 归 因 分 析 。 结 果 表 明 : 经 不 同 数 据 验 证 ,BEPS 模型 计算 结果 能 够 精确 反应 研究 区 


域 ET 及 其 组 分 的 分 布 情 况 和 变化 趋势 。 


1981—2018 年 研究 区 草地 .农田 和 和 森林 多 年 平均 ET 分 别 


为 278.22 mm .362.50 wa Oa mm, E,TfRET £ 4E E PELEAR, LARD A H 0.42 


mm*a1.0.63 mm:a 471.05 mm:a ^, ET 5 TEA A Z 


主要 受到 7 年 际 波动 的 影响 。 


田 7 和 ET 以 及 森林 的 ET 主要 受到 饱和 水 汽 压 差 (VPD ) 和 平均 气温 (TEMP) 变 化 的 控制 。 
森林 归 一 化 植被 指数 (Normalized difference vegetable index, NDVI) 4 © J/h 


NDVI AY Bh ERR, AAR 
X 键 词 : ABM; AA; 
文章 编号 : 


CEA A. 
蒸腾 ; BEPS 模 型 ， 半 干旱 区 


Z% BC (Evapotranspiration ,ET) 表 示 陆 面 水 分 转 
化 为 水 葵 汽 后 进入 大 气 的 过 程 … ,主要 由 蒸发 
(Evaporation , E ) fI12& f$ (Transpiration ,7 组 成 ,是 生 
态 系统 水 循环 中 的 重要 一 环 ,与 水 分 物理 传输 过 程 
和 能 量 转换 过 程 密切 相关 ”。 随 着 气候 变化 和 人 类 
城市 化 发 展 ,ET 发 生 了 明显 上 且 复 杂 的 变化 ”。 从 区 
域 尺 度量 化 蒸发 .蒸腾 和 蒸 散 情况 ,认识 环境 因素 
对 其 的 影响 机 制 ,有 助 于 科学 认识 区 域 水 文 循环 过 
程 ,促进 区 域 水 资源 管理 ,为 气候 变化 影响 生态 系 
统 水 文 循环 的 相关 人 研究 提供 科学 参考 “。 

半 干 旱地 区 生态 环境 较 脆弱 ,气候 变化 与 人 类 


活动 极 易 对 生态 系统 带 来 不 可 道 影响 ”“。 由 于 降水 
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间 分 布 格局 相似 ,与 五 相反 ,ET 年 际 波动 


综合 影响 因子 的 变化 和 忆 .T.ET 对 因子 的 敏感 性 ,研究 区 域 草地 和 和 农 


农田 和 
趋势 ,但 森林 环境 了 对 


少 ,水 资源 时 空 分 布 不 均匀 ,使 得 水 分 成 为 制约 半 
干旱 地 区 生态 系统 发 展 的 重要 因素 “"”。 蒸 散 作为 
半 干 旱地 区 水 循环 中 的 主要 组 成 之 一 ,其 变化 特征 
受到 植被 自身 特性 和 环境 条 件 的 影响 ,具有 显著 的 
区 域 异 质 性 和 季节 差异 性 ”"。 众 多 学 者 对 半 干 旱地 

区 葵 散 时 空格 局 和 变化 规律 进行 了 研究 。 代 鹏 超 
等 "基于 陆地 表面 能 量 平衡 算法 (Surface energy bal- 
ance algorithm for land, SEBAL) JF ZZ T Ti 10] jii MAR 
散 的 时 空 变化 ,认为 近 60 a 蒸 散 呈 波动 减 小 趋势 , 且 
空间 分 布 与 地 表 利 用 类 型 有 关 。 金 学 杰 等 "基于 
地 表 能 量 平衡 系统 (Surface energy balance system, 
SEBS) 模 拟 了 黑河 下 游 蒸 散 变 化 ,表明 不 同 地 表 利 
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与 气象 要 素 的 关系 ,认为 下 垫 面 供水 和 实际 水 汽 压 
变化 使 得 蒸 散 增 加 。 王 思 如 等 ”利用 生态 水 文 模 
型 ,分 析 了 科尔沁 沙 地 气象 因子 变化 和 蒸 散 的 关 
系 , 认 为 降水 的 年 际 变化 是 影响 蒸 散 的 主要 因子 。 
对 比分 析 ,发 现 我 国 半 干旱 地 区 蒜 散 空间 分 布 不 均 
5] ,变化 趋势 及 其 主 控 因 子 存 在 一 定 差异 FF RUP 
有 研究 将 蒸 散 拆 解 为 营 腾 和 莹 发 ,进一步 分 析 莹 散 
过 程 中 水 分 的 分 配 问 题 。 

气候 变化 与 人 类 活动 加 剧 的 背景 下 ,内 蒙古 半 
干旱 区 的 锡林郭勒 和 马兰 察 布地 区 暧 干 化 使 得 当 
地 水 分 供需 平衡 存在 不 确定 性 ”。 分 析 其 菩 散 时 
空 变 化 不 仅 能 帮助 揭示 水 文 循环 机 理 , 又 有 利于 评 
佑 区 域 水 分 气 缺 状况 ,为 当地 干旱 灾害 评估 与 水 资 
源 开 发 利用 提供 科学 依据 。 因 此 本 研究 选择 位 于 
半 干 旱 区 的 锡林郭勒 盟 和 乌 兰 察 布 市 作为 研究 区 ， 
分 析 和 量化 目标 区 域 蒸 发 .蒸腾 和 蒸 散 的 时 空 变化 
及 其 组 分 关系 ,对 比 不 同 气象 要 素 和 植被 生长 状况 
对 蒸发 .蒸腾 和 燕 散 变化 的 影响 ,揭示 半 干 旱 区 生 
态 系 统 燕 发 .蒸腾 和 蒸 散 对 气候 变化 的 响应 特征 及 
ERAT. 


1 材料 与 方法 


1.1 研究 区 概况 

本 文 研 究 区 域 为 内 蒙古 半 干 旱 区 的 锡林郭勒 
盟 与 马兰 察 布 市 ,地 理 位 置 位 于 40.16°~46.77°N， 
110.33°~119.91°E 之 间 ( 图 1)。 该 区 域 面 积 约 为 
28.56x10* km2? ,多 年 平均 总 降水 量 为 314.76 mm, & 
年 平均 气温 为 3.12 C ,属于 温带 大 陆 性 气候 ,地 表 
覆盖 类 型 主要 以 草地 和 农田 为 主 。 研 究 区 内 2 个 
代表 性 草地 生态 观测 站 "为 多 伦 站 和 锡林浩特 
站 ,经 纬度 分 别 为 42.53°N, 116.22°E Fil 43.53°N, 
116.67 灾 ,有 长 期 水 汽 通 量 观测 试验 ,用 于 本 研究 对 
比 验 证 。 

根据 地 表 植 被 覆盖 类 型 数据 统计 ,1981 一 2018 
年 研究 区 草地 、 农 田 和 森林 多 年 平均 分 别 占 
75.07% 11.67% 、5.94% ,植被 覆盖 分 布 面积 变化 较 
小 ,其 中 草地 面积 增长 1.50% ,农田 面积 减少 0.84%， 
森林 面积 共 减 少 1.56%。 因 此 本 人 研究 忽略 由 于 植被 
覆盖 类 型 变化 的 影响 ,主要 探究 气候 变化 和 植被 生 
长 状况 对 研究 区 域 水 循环 的 影响 。 


图 例 
= 国界 
— BF 
一 盟 市 边界 国 居民 用 地 
O 研究 区 — mm WU 


Em 草地 
gu kik 


E 农田 mm 森林 


注 :该 图 基于 国家 测绘 地 理 信 息 局 标准 地 图 服务 网 站 下 载 的 审 图 
号 为 GS(2019)3333 号 的 标准 地 图 制作 , 底 图 边界 无 修改 。 下 同 。 
图 1 研究 区 示意 图 
Fig. 1 Distribution map of study area 


1.2 模型 介绍 

本 文采 用 生态 系统 生产 力 模拟 (Boreal ecosys- 
tem productivity simulator, BEPS ) 模 型 模拟 区 域 的 蒸 
腾 蒸发 和 蒸 散 。 模 型 最 主要 的 特点 是 将 冠 层 叶片 
分 为 阳 叶 和 阴 叶 ,通过 进行 气孔 导 度 的 积分 对 叶片 
尺度 的 Farquhar 瞬 时 光化学 模型 进行 时 空 尺度 转 
换 ,模拟 计 算 逐 日 的 碳水 循环 通 量 ”“。 经 过 不 断 
发 展 和 完善 ,BEPS 被 证 实 能 够 准确 的 估算 生态 系统 
蒸 散 过 程 ,并 成 功 应 用 于 东亚 ”和 中 国 ” 等 地 。 
模型 的 ET 主要 计算 过 程 详 见 Chen M Liu E, 
1.3 数据 来 源 

模型 需要 的 气象 输入 数据 来 自 于 国家 气象 科 
学 数据 中 心 的 中 国 地 面 气候 资料 日 值 数据 集 (http:// 
data.cma.cn) ,包括 1981 一 2018 年 国家 气象 站 点 的 
气温 、 降 水量、 相对 湿度 .日 照 时 数 等 日 数据 。 该 数 
据 集 经 严格 的 质量 控制 ,包括 空间 和 时 间 上 的 一 致 
性 检查 及 相应 的 人 工 核 查 与 更 正 , 进 一 步 对 数据 进 
行 异 常 值 剔 除 。 研 究 区 域 地 势 较 平坦 ,因此 采用 
Kriging 方 法 插值 至 0.01° 并 根据 研究 区 域 边界 掩 膜 
裁 芍 , 绪 得 958x662 像 元 数 的 空间 日 数据 集 。 

1981 一 2018 年 的 叶 面 积 指 数 (Leaf area index, 
LAIT) 数 据 来 源 于 中 国 科学 院 地 理科 学 与 资源 研究 
所 刘 荣 高 团队 制作 的 每 8 d 的 1 km GLOBMAP LAI 
V3 产品 ,该 产品 经 验证 具有 较 高 的 精度 ”。 将 数据 
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韩 典 尾 等 ， 内 蒙古 半 干 旱 区 燕 散 估算 和 归 因 分 析 


重 采样 至 像 元 分 辨 率 为 0.01°, 并 掩 膜 至 与 气象 数据 
相同 大 小 ,统一 作为 模型 输入 数据 。 

地 表 利 用 类 型 数据 来 源 于 中 国 科 学 院 地 理科 
学 与 资源 研究 所 资源 环境 科学 数据 中 心 (http://www. 
resdc.cn/) 。 该 数据 集 以 Landsat TM/ETM 遥感 影像 
为 主要 数据 源 ,通过 人 工 目 视 解 译 生成 ,是 目前 我 
国 精度 最 高 的 土地 利用 遥感 监测 数据 产品 之 一 。 
土地 利用 类 型 主要 划分 为 农田 、 和 森林 草地、 水 域 、 
居民 用 地 和 裸 地 6 个 一 级 类 型 。 数 据 分 辨 率 为 
0.01° ,根据 研究 区 域 进行 裁剪 掩 膜 ,作为 模型 地 表 
徐 盖 类 型 输入 数据 。 

归 一 化 植被 指数 (Normalized difference vegeta- 
tion index, NDVI) 是 反映 植被 生长 状况 的 重要 参 
数 。 本 文选 用 1981 一 2015 年 GCIMMS3g 数 据 (https:// 
ecocast.arc.nasa.gov/data/pub/gimms ) 和 1998—2018 
年 中 国 科学 院 SPOT/VEGETATION NDVI 合 成 数据 
(http://www.resdc.cn)。 数 据 预 处 理 包括 对 原始 数据 
进行 格式 转换 .投影 转换 .数据 拼接 、 裁 切 等 ,并 将 
数据 重 采 样 为 0.01" 分 辩 率 的 影像 。 为 获得 连续 一 
致 的 NDVI, 使 用 最 大 值 合成 法 计算 NDVI 年 数据 *， 
用 于 区 域 蒸 散 及 组 分 变化 的 归 因 分 析 。 

站 点 验证 数据 采用 研究 区 内 多 伦 站 和 锡 林 浩 
特 站 的 开路 涡 度 相关 系统 观测 的 水 汽 通 量 数据 。 
利用 垂直 风速 脉动 的 平均 协 方差 和 水 汽 浓度 标量 
的 乘积 计算 30 min 平 均 潜 热 通 量 ,利用 EddyPro 软 
件 对 数据 进行 处 理 计算 ,包括 去 除 野 点 、 二 维 坐 标 
旋转 订正 等 。 将 每 日 潜 热 通 量 除 以 汽化 潜 热 (A~ 
2.45 MJ"kg') 计 算得 出 日 获 散 数据 。 对 于 缺失 时 间 
短 (<2 hh) 的 数据 采用 线性 内 插 法 进行 插 补 ,对 于 缺 
失 时 间 较 长 的 数据 采用 日 平均 法 对 其 进行 插 补 ”， 
通 量 数 据 的 能 量 平衡 比率 大 于 90% ,符合 能 量 不 闭 
合 程度 的 要 求 , 说 明 数 据 质量 良好 ” , 通 量 数据 的 具 
体 处 理 方法 介绍 详 见 Tian 等 的 文章 ””。 

1.4 敏感 性 分 析 

本 文采 用 由 MeCuen 提出 的 敏感 性 系数 , 即 
变量 与 单个 自 变 量 要 素 相 对 变化 量 之 比 , 反 映 各 自 
变量 要 素 对 因 变 量变 化 的 贡献 情况 ,具体 计算 方法 
见 毕 彦 杰 等 ”的 研究。 结合 敏感 性 系数 及 自 变 量 
本 身 变化 情况 ,综合 解析 和 量化 要 素 对 因 变 量 的 作 
用 。 本 文 主要 选择 平均 气温 (TEMP) 、 降 水 量 
(PRE) 饱和 水 汽 压 差 (VPD) 太阳 总 辐射 (RAD ) 和 
NDVI 5 个 要 素 分 析 其 变化 对 蒸发 . 藻 腾 和 蒸 散 变 化 


2 结果 与 分 析 


2.1 模型 验证 

图 2 给 出 了 BEPS 模型 模拟 的 研究 区 多 伦 站 和 
锡林浩特 站 2006—2008 年 月 蒸 散 与 相应 实测 值 对 
比 结果 。BEPS 模 型 月 ET 模拟 值 与 2 站 点 月 实测 值 
差距 较 小 ,RMSE 分 别 为 9.62 mm 和 10.14 mm (图 
2a~b)o BEPS 模型 可 以 解释 2 站 85% 的 蒸 散 变 化 


140 r (a) 多 伦 站 
120 上 RMSE( 模 拟 - 实 测 )-9.62 


一 一 ET 模拟 值 
一 6 一 ET 实测 值 


140 T (o) 锡林浩特 站 
120 上 RMSE( 模 拟 -实测 )=10.14 


— —— ET 模拟 值 
—e— ET 实测 值 


日 期 (年 -月 ) 
100 | (c) 实测 值 与 模拟 值 对 比 


y=0.86x 
R?-0.85 


月 ET 模拟 值 /mm 


20 40 60 80 100 
月 ET 实测 值 /mm 


注 :ET 为 蒸 散 。 下 同 。 
图 2 2006 一 2008 年 月 藻 散 模拟 值 和 实测 值 变 化 及 对 比 
Fig.2 Changes and comparison of simulated and measured 


values of monthly evapotranspiration from 2006 to 2008 
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X1 BEPS 模型 模拟 值 与 实测 值 数据 对 比 


Tab.1 Comparisons of simulated values from BEPS model and measured data 


纬度 /RN 经 度 /?E 年 份 ET 实测 值 /mm ET 模拟 值 /mm ”7 实测 值 /mm 7 模拟 值 /mm 7/ET 实 测 值 TET 模 拟 值 参考 文献 
42.05 116.28 2006 433.53 423.60 - - - = 30] 
43.55 116.67 2006 281.57 317.76 = = = - 30] 
44.13 116.33 2004 323.14 294.12 - - - - 31] 
2006 318.04 297.60 - - - - 31] 
43.53 116.67 2006 307.70 316.92 = = = - 32] 
2011 415.12 389.25 - - - - 32] 
43.63 116.71 2004 - - 155.25 158.64 - - 33] 
2006 = - 95.70 99.80 - = 33] 
43.60 116.70 2003—2005 = = = = 0.44 0.45 34] 
43.50 116.70 2003(5 一 9 月 ) = = - - 0.62 0.63 35] 
注 :ET 为 蒸 散 ;7 为 蒸腾 ;77ET 为 蒸腾 和 蒸 散 的 比值 。 下 同 。 
(图 2c)。 通 过 与 他 人 文献 结果 对 比 ( 表 1),BEPS 模 — 29-12.43 mm, 


型 模拟 了 时 平均 相对 误差 为 3.20% ,模拟 ET 时 平均 
相对 误差 为 6.70% ,模拟 TVET 时 平均 相对 误差 为 
1.90%, 这 表明 BEPS 模 型 在 站 点 尺度 模拟 准确 性 较 
好 23。 将 T ET 和 TET 空 间 分 布 情况 与 Niu 0" 
发 布 的 1981 一 2015 年 中 国 陆 地 生态 系统 蒸腾 蒸 散 
比 数据 集结 果 相 对 比 ,三 者 空间 分 布 情况 高 度 一 
致 。 综 上 ,BEPS 模 型 估算 结果 合理 ,可 以 用 于 研究 
区 ET 和 ET 的 模拟 计算 。 
22 蒸 散 及 其 组 分 的 时 空 特征 

根据 地 表 利 用 类 型 数据 ,提取 草地 农田 和 和 森 
林 像 元 ,并 统计 不 同 地 表 利 用 类 型 下 .7T 和 ET。 表 
2 给 出 了 研究 区 域 1981 一 2018 年 不 同 地 表 类 型 年 平 
均 耗 水 量 和 产 水 量 (PRE-ET) 的 总 体 情 况 。 草 地 的 
五 了 和 ET 耗 水 量 分 别 为 170.20 mm, 108.02 mm 和 
278.22 mm。 农 田 的 已 7 和 ET 耗 水 量 分 别 为 149.73 
mm 212.77 mm 和 362.50 mm. ZEB ET TI ET TÉ 
水 量 分 别 为 119.17 mm 189.64 mm 和 308.81 mm. 
结果 表明 ,农田 总 耗 水 量 最 大 ,森林 次 之 ,草地 最 
小 。 结 合 降水 量 看 ,和 森林 PRE-ET 最 大 ,为 53.52 
mm; 草地 次 之 ,为 31.32 mm; 农 田 PRE-ET 为 负 值 ， 


表 2 1981—2018 年 研究 区 不 同 地 表 类 型 年 平均 耗 水 量 
Tab. 2 Average annual water consumption of different 


surface types in the study area from 1981 to 2018 /mm 


地 表 类 型 ” PRE E 7 ET PRE-ET 
草地 309.54 170.20 108.00 27822 31.32 
农 350.07 149.73 212.77] 362.50 -1243 
森林 362.33 119.17 189.64 308.81 53.52 


注 :PRE 为 降水 量 ;E 为 蒸发 ;PRE-ET 为 产 水 量 。 


图 3 为 1981 一 2018 年 研究 区 ET 和 ET 的 年 际 
变化 。 结 果 表 明 研究 区 五 了 和 ET 多 年 平均 值 分 别 
为 167.13 mm 119.16 mm 和 286.35 mm, ZHI E 
著 上 升 趋势 (P<0.05) ,上升 速率 分 别 为 0.42 mm-a", 
0.63 mm-a'! fl 1.05 mm-a', 49 TAHE , E EREE 
较 平 稳 ,波动 幅度 更 小 ,因此 ET 年 尺度 波动 主要 受 
到 7 影响 ,变化 情况 与 了 更 一 致 。TVET 多 年 平均 值 
为 41.53% ,38 a 共 上 升 3.80% ,表明 蒸腾 在 蒸 散 中 的 
比例 时 增加 趋势 。 

图 4 为 研究 区 1981 一 2018 年 多 年 年 均 E、T、ET 
All TIET 的 空间 分 布 。E 多 年 平均 空间 分 布 整 体 呈 
西北 向 东南 递减 趋势 ,在 西北 部 植被 稀疏 区 五 较 高 ， 


350r. gp ..T y=1.05x+265.88, P<0.05 7 60 

300 L — REHE 
E 50 
Ri 250 上 5=0.06x440.32, P>0.05 Y ¢ " 
B V A wa \A : i 3 & 
S 200- Y B 
E y=0.42x+158.94, P<0.05, È 
Ej 150. CV si 

k eu A] — 9d 
100 f =0.63x+106.88, P«0.05 
1 ï 1 1 1 1 1 1 1 20 


1980 1985 1990 1995 2000 2005 2010 2015 2020 
年 份 


注 :五 为 蒸发 ;7 为 蒸腾 ;7T7ET 为 蒸腾 和 蒸 散 的 比值 。 下 同 。 
图 3 1981 一 2018 年 研究 区 蒸发 .蒸腾 . 蒸 散 蒸腾 和 
蒸 散 比 的 变化 趋势 


Fig.3 Trends of evaporation, transpiration, 


evapotranspiration and transpiration/evapotranspiration 


in the study area from 1981 to 2018 
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Fig.4 Spatial distributions of evaporation, transpiration, evapotranspiration and 


transpiration/evapotranspiration in the study area from 1981 to 2018 


Hf 35 160 mm, 但 在 东北 部 地 区 五 较 低 ,小 于 140 
mm。7 了 空间 分 布 和 EA 相反 ,在 西北 部 草地 稀 玖 区 7 
较 低 ,东北 部 和 东南 部 森林 农田 覆盖 区 较 高 。 在 
西北 部 地 区 ,7 小 于 80 mm ,在 东北 部 至 南部 的 森林 
农田 覆盖 区 ,7 大 于 120 mm。ET 和 TVET 空 间 分 布 
与 ?相似 ,草地 覆盖 区 ET 多 小 于 275 mm ,南部 和 东 
北部 农田 和 森林 覆盖 区 域 ET 上 且 大 于 300 mmo T/ET 
在 西北 部 植被 覆盖 度 较 低 区 小 于 0.40 ,南部 农田 履 
盖 区 和 东北 部 森林 覆盖 区 大 于 0.50。 

图 5 为 1981 一 2018 年 研究 区 TET 和 T/ET 的 
变化 趋势 空间 分 布 。E 在 全 区 域 多 呈 增 大 趋势 , 特 
别 在 东北 部 和 森林 和 中 南部 农田 区 呈 增 大 速率 且 大 
于 0.50 mm'a ,在 西北 部 和 南部 部 分 地 区 呈 减 小 趋 
势 , 但 减 小 趋势 不 显著 。 在 农田 覆盖 的 中 南部 和 南 


部 地 区 ,7 了 增 大 趋势 显著 (P<0.05) , 增 大 速率 大 于 
1.00 mma ,但 在 东北 部 森林 覆盖 区 ,了 呈 减 小 趋 
热 , 减 小 速率 最 快 为 -0.78 mma, A ERIT AE 
化 ,ET 在 全 区 域 多 呈 增 大 趋势 ,特别 在 北部 草地 覆 
六 区 和 南部 农田 覆盖 区 增 大 趋势 较 显 著 (P<0.05)， 
增 大 速率 大 于 1.50 mm'a3;ET 仅 在 西北 部 和 东北 部 
部 分 地 区 呈 减 小 趋势 ,但 不 显著 。TET 空 间 变 化 趋 
势 和 ?7 相似 ,从 整体 上 看 西部 地 区 多 呈 上 升 趋势 ,但 
在 东北 部 森林 覆盖 区 域 7ET 呈 下 降 趋势 , 减 小 速率 
最 快 为 -0.002.a'。 
2.3 蒸 散 及 其 组 分 的 归 因 分 析 

表 3 为 1981 一 2018 年 研究 区 主要 影响 因子 变化 
情况 。NDVI 区 域 平均 值 为 0.22 ,表现 为 增加 趋势 , 变 
化 率 为 0.0001.a- ,草地 呈 增 大 趋势 ,农田 和 森林 均 呈 
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图 5 1981 一 2018 年 研究 区 蒸发 .蒸腾 STL TPS RU ZR HA EC 2e HE ns [n] 
Fig.5 Spatial distributions of change rates for evaporation, transpiration, evapotranspiration 


and transpiration/evapotranspiration in the study area from 1981 to 2018 


表 3 1981 一 2018 年 研究 区 主要 影响 因子 的 变化 率 统计 


Tab.3 Change statistics of main influencing factors of the study area from 1981 to 2018 


地 表 类 型 NDVI/a" TEMP/C a" PRE/mm-:a' VPD/kPa-a* RAD/W m”. a” 
全 区 域 0.0001 0.05 -0.88 0.002 -0.11 
草地 0.0002 0.05 -0.92 0.002 -0.10 
农田 -0.0003 0.05 -0.32 0.002 -0.17 
森林 -0.0004 0.05 -1.00 0.002 -0.08 


注 :NDVI 为 归 一 化 植被 指数 ;TEMP 为 平均 气温 ;VPD 为 饱和 水 汽 压 差 ;RAD 为 地 表 总 辐射 。 下 同 。 


减 小 趋势 ,其 中 森林 减 小 趋势 最 大 ,为 -0.0004a'。 
TEMP 在 全 区 域 平均 为 3.12 "*C, 呈 增 大 趋势 ,变化 率 
为 0.05 °C+a',38 a 升 高 1.90 °C, PRE 在 全 区 域 多 年 
年 均 为 314.76 mm, 全 区 域内 呈 减 小 趋势 ,森林 减 小 
最 快 ,变化 率 为 -1.00 mm-a',38a FI f 38.00 
mm。VPD 在 全 区 域 多 年 年 均 为 0.37 kPa, 全 区 域 多 


时 增 大 趋势 ,平均 变化 率 为 0.002 kPa*a,38 a EFH 
0.09 kPa。RAD 全 区 域 呈 减 小 趋势 ,农田 减 小 速率 
最 快 ,变化 率 为 -0.17 W-*m^-a*,38a 下降 6.46 Wem’, 

1981—2018 年 各 要 素 变 化 对 区 域 五 变化 的 贡 
献 情况 如 表 4 所 示 。 和 森林 EE 变化 最 大 ,相对 变化 率 
J 15.90% ,考虑 五 对 各 要 素 变化 的 敏感 性 系数 及 各 
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要 素 相 对 变化 率 ,VPD 增 大 为 森林 EE 增加 的 主要 原 
,VPD 增 大 引起 E 增 大 12.22%。 草 地 和 农田 E 相 
对 变化 率 分 别 为 9.38% 和 10.90%。 与 森林 主导 因子 
相似 ,草地 和 农田 变化 的 主导 因子 也 为 VYPD。 综 
合 全 区 域 E 相 对 变化 率 为 9.55%,E 的 主导 因子 排序 
为 VPD>PRE>RAD>TEMP>NDVI。 

1981 一 2018 年 各 要 素 变 化 对 区 域 了 变化 的 贡献 
情况 如 表 5 所 示 。 农 田 7 相 对 变化 率 最 大 ,为 
24.80%。 综 合 7 对 各 因子 变化 的 敏感 性 和 各 因子 相 
对 变化 ,VPD 和 TEMP 增 大 为 农田 7 增 大 的 主要 原 
因 。 和 森林 7 了 变化 主要 受到 VPD 增 大 的 正 向 影响 和 
NDVI 减 小 的 负 向 影响 ,相对 变化 率 为 13.83%。 草 
地 了 相对 变化 率 为 19.34% ,主导 因子 和 农田 相似 。 
综合 全 区 域 ,T 相 对 变化 率 为 20.09% ,主导 因子 依次 
为 VPD>TEMP>NDVI>RAD>PRE。 

1981 一 2018 年 各 要 素 变化 对 区 域 ET 变 化 的 贡 
献 情 况 如 表 6 所 示 。 农 田 ET 相 对 变化 率 最 大 ,为 


19.10%。 各 要 素 对 农田 ET 相对 贡献 大 小 依次 为 
VPD>TEMP>RAD>NDVI>PRE, VPD 和 TEMP 增 大 
为 农田 了 T 增 大 的 主要 原因 。 草 地 和 森林 相对 变化 
率 分 别 为 13.25% 和 14.64% ,各 要 素 对 草地 和 森林 
ET 相对 贡献 依次 为 VPD>TEMP>NDVI>RAD> 
PRE。 综 合 全 区 域 了 并 相对 变化 率 为 13.93% ,主导 
子 排序 为 VPD>TEMP>RAD>PRE>NDVI。 


3 讨论 


人 研究 通过 站 点 ET 实测 数据 和 前 人 模拟 结果 相 
对 比 验 证 发 现 ,T、ET 和 TET 售 算 相 对 误差 最 小 分 
别 为 2.10% 、2.30% 和 1.60%。Li 等 中 研究 结果 表 
HH ,TET 和 LAI 水 平 相关 , 当 LAI 较 小 时 ,7/ET 稳 定 
在 0.33+0.50 , 而 在 植被 生长 状况 良好 区 域 ,TVET 大 
于 0.40。 本 研究 区 西北 部 LAI 数 值 较 小 ,变化 范围 
为 1.30~2.00 m? - m?, 4 KREI TIET X 0.38 , ifii Zi 


表 4 1981 一 2018 年 研究 区 影响 因子 对 蒸发 变化 的 贡献 分 析 


Tab.4 Contribution analysis of influencing regulators to evaporation dynamics in the study area from 1981 to 2018 


"e 相对 贡献 /% 
地 表 类 型 EE 相对 变化 率 /9% R 主导 因子 排序 
NDVI TEMP PRE VPD RAD 
全 区 域 -0.08 1.11 1.94 732 | -132 9.55 0.94 VPD>PRE>RAD>TEMP>NDVI 
草地 -0.17 1.09 1.50 7.38 | -072 9.38 0.96 VPD>PRE>TEMP>RAD>NDVI 
农田 1.45 1.11 1.02 8.54 -1.84 10.90 0.94 VPD>RAD>NDVI>TEMP>PRE 
森林 1.46 2.25 1.19 12.22 -2.46 15.90 0.92 VPD>RAD>TEMP>NDVI>PRE 


注 :5 个 自 变量 要 素 相 对 贡献 之 和 等 于 因 变 量 燕 发 相对 变化 乘 以 尼 ; 尼 为 整体 拟 合 优 度 。 下 同 。 


表 5 1981 一 2018 年 研究 区 影响 因子 对 蒸腾 变化 的 贡献 分 析 


Tab.S Contribution analysis of influencing regulators to transpiration dynamics in the study area from 1981 to 2018 


相对 贡献 /% 


地 表 类 型 7 相对 变化 率 /% R 主导 因子 排序 
NDVI TEMP PRE VPD 
全 区 域 3.42 5.41 -1.12 1150 -2.37 20.09 0.84 VPD>TEMP>NDVI>RAD>PRE 
草地 4.67 5.35 -1.63 1013 -1.74 19.34 0.86 VPD>TEMP>NDVI>RAD>PRE 
农田 -1.31 8.37 -1.12 1649 -1.63 24.80 0.84 VPD>TEMP>RAD>NDVI>PRE 
森林 -3.49 3.27 -0.63 1304  -144 13.83 0.85 VPD>NDVI>TEMP>RAD>PRE 


426 1981 一 2018 年 研究 区 影响 因子 对 蒸 散 变化 的 贡献 分 析 


Tab. 6 Contribution analysis of influencing regulators to evapotranspiration dynamics in the study area from 1981 to 2018 


i 相对 贡献 /% " . : 
地 表 类 型 ET 相对 变化 率 /% R 主导 因子 排序 
NDVI TEMP PRE VPD RAD 
全 区 域 0.79 4.96 -1.01 10.84 -2.84 13.93 0.91 VPD>TEMP>RAD>PRE>NDVI 
草地 2.44 2.88 -0.33 8.65 -1.49 13.25 0.91 VPD>TEMP>NDVI>RAD>PRE 
农田 -1.54 6.65 -0.65 15.34 -2.07 19.10 0.92 VPD>TEMP>RAD>NDVI>PRE 
森林 -1.32 3.27 —0.63 13.04 -1.14 14.64 0.90 VPD>TEMP>NDVI>RAD>PRE 
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部 森林 和 南部 农田 覆盖 区 植被 生长 状况 较 好 ,7T/ET 
则 稳定 在 0.50 左 右 ,与 Li 等 研究 结果 相符 。1981 一 
2018 年 BEPS 模拟 多 年 平均 ET 结果 介 于 278.22~ 
362.50 mm; 牛 忠 恩 等 2 基于 PT-JPL 模型 对 中 国 ET 
分 布 情况 进行 估算 ,内 蒙古 地 区 多 年 平均 ET 介 于 
300-400 mm; 王 思 如 等 所 基于 WaVEM 模型 对 内 蒙 
古 科 尔 沁 地 区 ET 变化 进行 研究 ,发 现 ET 介 于 200~ 
450 mm, 以 上 均 与 本 人 研究 模拟 结果 相符 。 因 此 通过 
站 点 实测 数据 及 研究 结果 的 比 对 证 明 BEPS 模 型 结 
果 合 理 , 可 以 较 准 确 模拟 研究 区 ET 及 其 组 分 。 

半 王 旱 区 不 同 的 生态 条 件 会 影响 茸 散 对 气候 
条 件 的 响应 吕 。 根 据 对 瑟 . 了 和 ET 空间 分 布 结果 UT 
究 区 西北 部 植被 稀 玻 ,年 均 NDVI 小 于 0.15 ,但 较 高 
的 气温 及 较 大 的 VPD 使 得 蒸发 速率 较 快 。 研 究 区 
东部 和 南部 森林 和 农田 环境 植被 覆盖 度 较 高 ,年 均 
NDVI 大 于 0.27, 下 垫 面 粗糙 度 更 高 ,植被 的 蒸腾 作 
用 成 为 地 表 和 大 气 间 水 分 交换 的 主要 方式 ,因此 在 
东北 部 森林 和 南部 农田 覆盖 区 域 T 和 7/ET 均 高 于 
其 他 区 域 。 根 据 Jiang 等 ”研究 , 半 干 旱 区 7 是 主导 
ET 变化 的 主要 驱动 组 分 ,结合 图 3 和 图 4,E 的 高 值 
区 主要 分 布 在 西北 部 植被 稀 玖 区 ,年 际 间 波 动 主要 
受到 气候 条 件 支 配 ,而 7 高 值 区 主要 分 布 在 东北 部 
森林 和 南部 农田 覆盖 区 ,年 际 变化 除 受到 气候 条 件 
影响 外 还 受到 植被 生长 状况 以 及 人 为 因素 (例如 农 
田 人 工 灌溉 等 ) 影 响 ,因此 五 年 际 波动 较 小 ,变化 趋 
势 较 平 稳 , 而 ET 和 7 年 际 波动 变化 一 致 性 较 高 。 

mVET 表 征 植被 燕 腾 对 生态 系统 蒸 散 的 贡献 率 ， 
是 准确 量化 生态 系统 水 分 利用 效率 的 关键 参数 ”。 
TEMP, VPD 和 NDVI 升 高 和 地 表 植 被 生长 会 影响 7/ 
ET 。TEMP 增 大 ,有 利于 延长 植被 生长 季 ;VPD 作 
为 大 气 水 分 气 缺 状况 的 表征 因子 ,主要 通过 影响 裸 
地 和 植被 表面 的 水 分 梯度 ,调节 叶片 气孔 导 度 影响 
植被 燕 腾 量 '…。 半 干旱 地 区 升温 效应 增 大 了 水 分 
梯度 和 叶片 气孔 导 度 ,加 速 叶 片 和 外 界 水 汽 交换 ， 
对 WET 起 促进 作用 。NDVI 增 大 表明 研究 区 域 植被 
长 势 逐 渐 增 加 , 冠 层 截获 的 辐射 增多 ,导致 蒸腾 耗 
水 增加 ,同时 长 势 良好 的 树木 和 作物 水 分 利用 效率 
更 高 , 燕 腾 作用 随 之 增强 。 结 合 E 和 7T 归 因 分 析 
结果 ,人 研究 区 草地 NDVI 多 年 呈 增 大 趋势 ,对 草地 
为 负 贡 献 ,对 7 为 正 贡 献 , 有 利于 增 大 草地 生态 系统 


燕 散 中 蔡 腾 的 比例 ,提高 了 草地 的 水 分 利用 效率 。 
而 森林 和 农田 NDVI 多 年 呈 下 降 趋 势 ,特别 在 森林 
敌 盖 区 NDVI 下 降 速 率 最 快 ,因此 农田 和 和 森林 部 分 
区 域 7/ET 呈 减 小 趋势 ,在 ET 中 的 占 比 逐 渐 增 加 。 

根据 前 人 研究 ,ET 和 气候 条 件 地表 和 覆盖 类 型 
有 密切 关系 。 如 表 2, 研 究 区 域 单位 面积 草地 万 
大 于 农田 和 森林 ,但 7 小 于 农田 和 森林 ,这 主要 由 于 
草地 植被 叶 面 积 较 小 ,地 表土 壤 裸 露 较 多 ,到 达 地 
表 的 光照 辐射 能 量 使 得 地 表 燕 发 旺盛 ,而 森林 和 农 
田 植被 叶片 面积 较 大 ,更 高 的 冠 层 高 度 和 植被 层 厚 
度 有 利于 吸收 更 多 光照 辐射 能 量 ,导致 下 渗 到 土壤 
表层 的 水 分 不 易 蒸发 ,进而 被 植被 发 达 的 根系 吸 
收 ,为 蒸腾 作用 提供 更 多 水 分 来 源 。 另 一 方面 , 农 
田 多 分 布 在 研究 区 域 南 部 ,平均 年 PRE 为 350.07 
mm, 但 PRE-ET 仅 为 -12.43 mm ,说 明 农 牧 资源 的 不 
合理 分 配 ,可 能 导致 未 来 半 干 旱 区 域 更 容易 发 生 圭 
壤 和 干旱 ,进而 造成 该 区 域 水 土 流失 及 土壤 退化 等 环 
境 问 题 。 


4 结论 


Exi 


本 文 以 位 于 内 蒙古 半 干 旱 区 的 锡林郭勒 盟 和 
乌 兰 察 布 市 为 例 ,结合 遥感 .气象 资料 和 站 点 观测 
资料 等 ,采用 BEPS 模型 ,分 析 和 量化 以 农田 .草地 
和 森林 为 主 的 半 干 旱 区 一 了 和 ET 的 时 空 变化 ,揭示 
半 干 旱 草 地 生态 系统 .7T 和 ET 对 气候 变化 的 响应 
特征 及 主 控 因子 ,得 出 主要 结论 如 下 : 
经 过 站 点 实测 数据 对 比 ,BEPS 模 型 计算 结果 能 
够 精确 反应 研究 区 ET 及 其 组 分 的 分 布 情况 和 变化 
趋势 。 研 究 区 域 草 地 农田 和 和 森 林 年 均 ET 分 别 为 
278.22 mm ,362.50 mm 和 308.81 mm, E,T fI ET £ 
年 旦 显著 上 升 趋势 ,上 升 速 率 分 别 为 0.42 mm-a, 
0.63 mm .a ! fll 1.05 mm-a/5,1981—2018 4E H E FH 
15.96 mm,23.94 mm 和 39.9 mm, ET 5j T E 4 [X bk 
内 空间 分 布 格局 相似 ,与 相反 ,ET 年 际 波动 主要 
受到 7 年 际 波动 的 影响 。 综 合影 响 因 子 的 变化 和 
TET 对 因子 的 敏感 性 ,人 研究 区 域 草地 和 农田 7 和 
ET 以 及 和 森林 的 ET 主要 受到 VPD 和 TEMP 变 化 的 控 
制 。 农 田 和 森林 NDVI 都 呈 减 小 趋势 ,但 森林 7 对 
NDVI 变 化 更 加 敏感 ,因此 负 贡 献 更 大 。 
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Abstract: Evapotranspiration (ET) is the most important flux term in the ecosystem water budget because it is 
related to water and heat exchange between ecosystems and the atmosphere. Quantifying the magnitude and 
dynamics of ET and its components include evaporation (E) and transpiration (7) at regional scales and deriving a 
mechanistic understanding of the underline regulations is essential for managing water resources under the 
changing climate. On the basis of a mechanistic ecosystem model (BEPS), we estimated the magnitudes of ET 
and its components in the semi-arid region (40.16?—46.77?N and 110.33?- 119.91? E) of Inner Mongolia, China 
from 1981 to 2018. Compared with different data, the BEPS could provide accurate estimates for changes in ET 
and its components by validating with ground measurements. The total average annual ET of grassland, cropland, 
and forests in the study area from 1981 to 2018 was 278.22 mm, 362.50 mm, and 308.81 mm, respectively. Over 
the study period, E, T, and ET increased at a rate of 0.42 mm-a ', 0.63 mm-a ', and 1.05 mm-a ', respectively. 
Across the study region, we found that the spatial distribution patterns of T was opposite to that of E and similar 
to that of ET. From the interannual variation standpoint, the interannual fluctuations in ET were primarily 
influenced by those in T. Considering the changes in influencing factors and the sensitivities of the three fluxes, 
we concluded that T and ET of the grasslands and croplands and ET of forests were mainly controlled by vapor 
pressure deficit and air temperature. Although both cropland and forest normalized difference vegetation index 
(NDVI) showed decreasing trends, forest environmental T was more sensitive to changes in NDVI, thus receiving 
more negative influence than cropland. 


Key words: evapotranspiration; evaporation; transpiration; BEPS model; semi-arid area 


